Sound focusing by gradient index sonic lenses 
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Gradient index sonic lenses based on two-dimensional sonic crystals are here designed, fabricated 
and characterized. The index-gradient is achieved in these type of flat lenses by a gradual modi- 
fication of the sonic crystal filling fraction along the direction perpendicular to the lens axis. The 
focusing performance is well described by an analytical model based on ray theory as well as by 
numerical simulations based on the multiple-scattering theory. 

PACS numbers; 43.20.Fn;43.58.Ls;43.20.Dk 



A two-dimensional (2D) sonic crystal (SC) is just a 
periodic distribution of solid cylinders in air with their 
axis parallel aligned. The properties of these structures 
have been already analyzed in a broad range of frequen- 
cies. Thus, for frequencies of the order of the lattice 
separation they present a band of frequencies (bandgap) 
in which the sound propagation is forbidden because of 
Bragg reflection [1]. However, in the range of low fre- 
quencies (homogenization limit) they behave like homo- 
geneous media whose effective acoustic parameters, dy- 
namical mass density and bulk modulus, basically de- 
pend on the lattice filling fraction[3-0|- 

The homogenization properties of SCs have been em- 
ployed to design refractive devices like, for example, 
acoustic lenses whose focusingproperties are based on 
their external curved surfaces [5|-l7| or Fabry- Perot like 
acoustic interferometers 0. More recently, gradient in- 
dex (GRIN)acoustic lenses based on SCs have been 
proposedjoifiol. Like their optical counterparts, the pro- 
posed GRIN SC lenses have flat surfaces and their index 
gradient is obtained by changing the SC filling fraction, 
which is directly related to the local index of refraction. 
An advantage of GRIN lenses in comparison with curved 
lenses is their easy fabrication. The bending of sound 
waves obtained by GRIN acoustics structures have been 
also proposed to produce acoustic mirages in a labfllj. 

In this letter 2D GRIN SC lenses have been fabricated 
and characterized by an specifically designed experimen- 
tal set up. Each lens is made of a rectangular cluster 
of aluminum cylinders distributed in an hexagonal lat- 
tice. Their focusing properties are here comprehensively 
studied either experimentally as well as theoretically. We 
demonstrate that their focusing behavior is a true GRIN 
effect that is fairly well reproduced by simulations based 
in the multiple scattering theory and it is also supported 
by an analytical model based on ray theory. To the 
best of our knowledge, no experimental demonstration 
of GRIN SC lenses had been presented up to date. 

The acoustical axis of GRIN lenses considered here is 
defined along the x— axis and, therefore, the refractive 
index n{y) changes along the y— axis. The local varia- 
tion of n is obtained by changing the cylinders radii as 
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FIG. 1: Schematic view of a GRIN lens made of nine columns 
of metal rods. Circles define the rods sections, which are 
position dependent. The sound impinging the lens from the 
left bends towards the x axis and focuses at position xf (red 
dot). The gray rectangle defines the area where the data are 
acquired. 



explained in Ref. Q. A GRIN lens and the region of 
data acquisition are schematically shown in Fig. [T] 

The index profile n(y) is designed with the goal of hav- 
ing low aberration of the focal spot. We choose a hyper- 
bolic secant profile because of its demonstrated property 
of focusing sound inside large enough GRIN structures 
without aberration lldl: 



n{y) = no sech{ay), 



where a is a constant defined as 
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where h is the half-height of the lens, uq is the refractive 
index on the a;— axis [y —O) and nu that at the lens edges 
(2/ = ±h). 

The differential equation satisfied by the energy rays 
can be obtained 12] and by solving it, the position of the 
focal spot outside the lens, x^, is obtained 
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TABLE I: Effective refractive index profile rief / (y) of the 2D 
GRIN SC lenses studied here. 
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where d is the lens thickness. 

A 2D GRIN SC can be considered as a homogenized 
medium made of subwavelength discrete units. The units 
are metal cylinders that are distributed in rows of equal 
radii. A given row £ is parallel to the x— axis and it is in- 
dependently homogeneous and isotropic. The refractive 
index at row with coordinate yi is adjusted by taking 
into account that, in the low frequency limit, a cluster of 
rigid cylinders embedded in a fluid or gas background be- 
haves like a homogeneous medium whose effective refrac- 
tive index Uef / mainly depends on the fraction of volume 
occupied by the cylinder in the SC or filling fraction0,|3]. 

The 2D GRIN SC lenses studied are designed with a 
profile n,{y) as closer as possible to that in Eq. ([T|) with 
no = 1.337 and = 1. Each lens is made of a rect- 
angular cluster of aluminum cylinders distributed in an 
hexagonal lattice with parameter a — 2cm (see Fig. [1]). 
The radiii Ri of cylinders at row ye is obtained by impos- 
ing the condition Ueffif) = n{yt), where / = ^{Rg/a)'^ 
is the filling fraction of the hexagonal lattice. Table U 
reports the values of Ri obtained. Note that values cor- 
responding to the lower half of the SC lens are the only 
ones reported since the lens is symmetric with respect to 
the cc— axis. Odd layers have 11 cylinders while even lay- 
ers have 10 cylinders. Column n{yi) in Table report the 
exact hyperbolic secant profile while those under Uef / {yi) 
are those effectively achieved in samples. 

Let us point out that the chosen lattice parameter 
determines an upper frequency limit for the homoge- 
nization. Since the homogenization condition is accom- 
plished while A > 4a[il, where a =2 cm, for frequencies 
V > 'i.2kHz the homogenization condition is broken. 

The lenses have been characterized in a 2D anechoic 
chamber specifically designed for the parameters of the 
lenses. The chamber consists of two parallel metal sheets 
separated a distance H = 5cm. For the frequencies an- 
alyzed here, from 3.5 kHz up to 4.5 kHz, there are two 



modes propagating trough the chamber, but the second 
one has a much larger wavelength and is considered as 
an offset. The sound waves are generated with a UDE 
AC-150 column loudspeaker separated 2 meters from the 
chamber in order to get approximately a plane wavefront 
at the entrance. The input signal is a Gaussian pulse 
centered at the frequencies in between the range under 
study. A total of eleven frequencies are comprehensively 
analyzed; i.e., from 3.5 kHz till 4.5 kHz in steps of 0.1 
kHz. Measurements made without lens confirm that the 
propagating field inside the chamber has practically a 
plane wavefront. We use a B&K 4958 microphone to 
record the transmitted sound behind the lens and use the 
PCI Acquisition Card NI-5105 for processing the data in 
a computer. Data are acquired in a square area of 24x24 
cm^, which is represented in Fig. [T] The acquisition 
points are separated by equal distance along the x— and 
?/— axis; i.e.. Ax = Ay ~ 2cm. 

Seven samples have been experimentally characterized. 
The lenses differ in the number TV of layers defining 
their thickness along the a;— axis; from 4 to 10. The 
corresponding pressure maps have been measured and 
compared with those obtained by numerical simulations 
based on the multiple scattering theory. Besides, the 
focal spots Xj have been also located and their values 
are compared with those predicted by Eq. Mea- 
surements and simulations are performed in the range of 
frequencies explored every 0.1 kHz. 

Numerical simulations are performed by using a 2D 
multiple scattering algorithm as described in Refs. 0, H, 
where the reader is addressed for details. In brief, a 2D 
sound wave with a plane wavefront impinges the lens 
from the left and the pressure map is calculated in the 
opposite side. As a typical example of results obtained. 
Figs, mja) andlDJb) show the sound amplification maps 
for a 9 layer thick GRIN lens for the frequency of 4.5 kHz; 
their corresponding wavelength, A « 3.8 a, roughly cor- 
responds to the cutoff defining the homogenization limit. 
The focusing effects, focal spot and diffraction lobes, are 
clearly seen in both maps, measured [see Fig. [DJa)] and 
calculated [see^b)]. A better comparison between data 
and simulations is given in Figs. |3l^a) and[3l^b) that show 
the longitudinal and transversal profiles passing through 
the focal spot obtained from the previous pressure maps. 
Note the good agreement exists between data and sim- 
ulations. A similar agreement is also obtained for the 
other frequencies analyzed. 

In spite of the agreement reported above between the- 
ory and experiment a comparison with the position of 
the focal spot x/ predicted by the ray theory of Eqs. (1) 
to (3) has been also performed in order to support that 
sound focusing is a truly refractive effect produced by 
the bending of waves inside the GRIN lens. The param- 
eters used in Eq. dSj are determined as follows. The lens 
half-height h is obtained from 2/i = Nya where Ny is the 
maximum number of cylinders along the y— axis, =11. 



3 




1.5 



0.5 

I. 



1.5 



0.5 



10 15 20 
X[cm] 



[t] 



FIG. 2; (a) Sound amplification map generated by 4.5 kHz 
sound waves impinging a 9 layers thick 2D GRIN SC lens, (b) 
The corresponding map obtained by using a multiple scatter- 
ing algorithm. Maps are obtained in an area 24x24 cm^ near 
the lens (see Fig. [1]). The asterisks mark the focal spot Xf. 
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FIG. 3: (a) Longitudinal and (b) transversal cuts extracted 
from the maps shown in Fig. [51 Simulations based on multiple 
scattering theory (blue curves) and experimental data (black 
dots) are depicted. 




_ 1 


1 


- -•- - Exp. 


4.5 kHz (b) . 






- -□- - MST 








Ray Model 


























1 








4 


5 


6 7 8 
Lens thickness (N) 


9 10 



FIG. 4: Focal spot, a;/, as a function of A'^, the number of lay- 
ers defining the lens thickness, for frequencies: (a) 3.5kHz and 
(b) 4.5kHz. Black dots and hollow squares represent exper- 
imental data are results obtained from a multiple scattering 
theory (MST), respectively. The dashed lines between sym- 
bols are guides for the eye. Red lines represent predictions 
obtained by a ray model (see text). 



where N is the number of cylinders' layers employed in 
building the lens and k is an adjustable parameter used 
to fit Xf obtained from Eq. ([3]) to the data. The values 
obtained are k = -0.41cm for 3.5 kHz, and k — -0.28cm 
for 4.5 kHz, which means that the lens surface is nearer 
to the lens units than that usually employed in solid state 
physics 13], where the surface position is taken at a dis- 
tance above the atoms equal to a half of the lattice pa- 
rameter; i.e.. k —0 in our model. 

Figure H] reports results obtained for Xf at the lower 
and higher frequencies analyzed: (a) 3.5kHz and (b) 
4.5kHz. 

Data and multiple scattering simulations show a bet- 
ter agreement with the ray model when: i) the frequency 
operation is near the homogenization limit and ii) for the 
thicker lens. These conclusions are physically intuitive. 
On the one hand, for lower frequencies (larger wave- 
lengths) the homogenization takes place over a larger 
number of unit cells in the lattice and, as a consequence, 
the locally designed refractive index profile n{y) could be 
partially destroyed. On the second hand, thicker lenses 
means that sound travels longer distance inside the lens 
and consequently the role of surfaces and diffraction ef- 
fects become negligible. 

In summary, we have demonstrated the focusing prop- 
erties of 2D GRIN sonic lenses made of aluminum rods. 
Their broadband performance has been confirmed by nu- 
merical simulations based on a multiple scattering algo- 
rithm and by an analytical model based on ray theory. 
It can be concluded that GRIN sonic lenses are a new 
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type of refractive acoustical devices that are feasible and 
reliable for possible applications not only in the sonic 
range, as has been demonstrated here, but also for ultra- 
sonics and even for acoustic surface waves as has been 
previously proposed. 

Work partially supported by USA Office of Naval 
Research (Grant N000140910554) and by the Spanish 
MICIIN under contracts TEC2007-67239 and CSD2008- 
66. 



* |jsdehesa@upvnet.upv.es| 
[1] J.V. Sanchez-Perez et al, Phys. Rev. Lett., 80, 5325 
(1998). 

[2] A. A. Krokhin and J. Arriaga and L.N. Gumen, Phys. 

Rev. Lett. 91, 264302(2003). 
[3] D. Torrent, A. Hakansson, F. Cervera and J. Sanchez- 



Dehesa, Phys. Rev. Lett. 96, 204302(2006). 
[4] D. Torrent and J. Sanchez-Dehesa, Phys. Rev. B 74, 
224305 (2006). 

[5] E. Meyer and E. G. Neumann, Physical and Applied 
Acoustics, (Academic Press, New York, 1972). 

[6] F. Cervera et al, Phys. Rev. Lett., 88, 023902 (2002). 

[7] B.C. Gupta and Z. Ye, Phys. Rev. E 67, 036603 (2003). 

[8] L. Sanchis, A. Hakansson, F. Cervera, and J. Sanchez- 
Dehesa, Phys. Rev. B67, 035422 (2003). 

[9] D. Torrent and J. Sanchez-Dehesa, New J. Phys. 9, 323 
(2007). 

[10] S.-C. S. Lin, T.J. Huang, J.H. Sun and T.T. Wu, Phys. 

Rev. B 79, 094302 (2009). 
[11] S.-C. S. Lin and T.J. Huang, J. Appl. Phys., 106, 53529 

(2009). 

[12] C. Gomez-Reino, M. V. Perez, and C. Bao, Gradient- 
index Optics: Fundamentals and Applications, (Springer, 
Berlin,2002). 

[13] K.W.-K. Shung and Y.C. Tsai, Phys. Rev. B 48,11265 
(1993). 



